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ABSTRACT 



Recently, a smart antenna, i.e., a blind adaptive antenna 
array, has attracted much attention to suppress multiple 
access interference (MAI) and multipath signals to improve 
the capacity of a code division multiple access (CDMA) 
wireless communications system. Most of the smart antenna 
algorithms either need matrix computation or complex cal- 
culations. In this invented direction-of-arrival (DOA) track- 
ing algorithm only 2M complex multiplications are required 
per snapshot. Further more the DOA tracking algorithm is 
not sensitive to the mixer phase distortions. Hence the 
system complexity is reduced since a separate phase cali- 
bration is not required for the DOA tracking algorithm. Also 
the equivalent DOA tracking error due to the mixer phase 
distortions is derived. Simulation results show that the DOA 
tracking algorithm in the present invention works effectively 
under additive white Gaussian (AWGN) as well as Rayleigh 
fading environments. Also the system employing the 
invented smart antenna has the same bit error rate (BER) 
performance as those of other existing smart antennas, 
however, has less computation burdens and complexity. 

48 Claims, 9 Drawing Sheets 
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DIRECTION OF ARRIVAL ANGLE 
TRACKING ALGORITHM FOR SMART 
ANTENNAS 

BACKGROUND OF THE INVENTION 5 

1. Field of the Invention 

The present invention relates to wireless telecommunica- 
tions in general and more particularly to a method and 
system for determining a direction of arrival angle with a 
smart antenna in a wireless telecommunications system. 10 

2. Related Art 

A smart antenna, such as a blind adaptive antenna array, 
is a method and system used to enhance a desired user 
signal, to suppress multiple access interfering signals, and in 15 
general to improve the performance of a wireless commu- 
nications system, including CDMA systems such as third 
generation (3G) wireless communication systems, particu- 
larly those specified in the wide band CDMA(W-CDMA) or 
the CDMA2000 standards. In such systems, smart antenna 2 o 
technology may optionally be employed at a base station. A 
smart antenna can suppress interfering signals of different 
DOAs (Directions of Arrival) by using spatial diversity. 
Smart antenna technologies attract much attention these 
days as they support more users with a high quality of 25 
service and high data rates. Efficient smart antenna schemes 
have appeared recently. It is desirable that a smart antenna 
be effectively easy to implement. 

The DO A tracking capability, convergence, and compu- 
tation complexity are very important issues in the design of 30 
the smart antennas. Conventional methods for determining 
the direction of arrival of an incoming signal in a smart 
antenna system require matrix computation or complex 
computation, which are the costly parts of the smart antenna 
tracking algorithms. Specifically, known algorithms attempt 35 
to calculate the weight vectors directly by minimizing or 
maximizing some cost criteria such as SINR (signal to 
interference+noise ratio), BER (bit error rate) or MSE (mean 
square error). These methods require determining a separate 
weight vector for each antenna element; e.g., if there are six 40 
antenna elements, six unknown variables must be deter- 
mined. Another drawback associated with known methods is 
that practical embodiments of these methods require cali- 
bration of mixer phase distortions. 

SUMMARY OF THE INVENTION 45 

The aforementioned needs are addressed to a great extent 
by the present invention, which provides a method and 
system for performing an initial DOA estimate and fine- 
tuning the estimate and for tracking the DOA when there is 50 
relative motion between the transmitter and receiver. The 
present invention utilizes the notion that the maximum 
antenna gain can be achieved when the weight vector for the 
antenna array is matched to the array response vector. In 
highly preferred embodiments, rather than attempting to 55 
calculate the weight vectors directly as is done in conven- 
tional algorithms, the present invention estimates a DOA 
angle and calculates the weight vectors as a function of the 
estimated DOA angle. The estimated DOA angle is then 
adaptively updated to optimize correlation output, e.g., to 60 
achieve the maximum output power or maximum signal-to- 
interference-plus-noise-ratio (SINR) for the desired user. 
This method may be used during both the initial search for 
the DOA, in which a coarse adjustment is made to the 
estimated DOA angle, and for tracking the DOA as relative 65 
movement between the transmitter and the receiver occurs, 
in which a fine adjustment is made. The proposed algorithm 
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of the present invention is not sensitive to mixer phase 
distortions. Thus, a separate phase calibration is not required 
in the system. The proposed algorithm also has lower 
computational burdens, as compared to conventional 
algorithms, for calculating weight vectors for smart antenna 
arrays, thereby lowering the complexity and cost of a 
system. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The aforementioned advantages and features of the 
present invention will be more readily understood with 
reference to the following detailed description and the 
accompanying drawings in which: 

FIG. 1 is a block diagram of a complex pseudo-noise (PN) 
spread receiver for a CDMA communication system con- 
figured in accordance with one embodiment of the present 
invention. 

FIG. 2 is a flow chart of a proposed adaptive DOA 
tracking algorithm in accordance with one embodiment,of 
the present invention. 

FIG. 3 is a beam pattern for an M-6 elements antenna 
array system and DOA=0°, configured in accordance with 
one embodiment of the present invention. 

FIG. 4 is a beam pattern for an M=2 elements antenna 
array system and DOA=60°, configured in accordance with 
one embodiment of the present invention. 

FIG. 5 is a graph of the probabilities of correct direction 
detection at the initial search versus the bit-energy-to-noise 
power spectrum density E^/Nq in dB for antenna arrays of 2, 
4, and 6 elements, configured in accordance with one 
embodiment of the present invention. 

FIG. 6 shows simulation DOA tracking results versus 
iteration with averaging over 100 runs under an additive 
white Gaussian noise (AWGN) channel. 

FIG. 7 shows simulation DOA tracking results versus 
iteration with averaging over 100 runs under a Rayleigh 
fading environment for one embodiment of the present 
invention. 

FIG. 8 shows simulation BER results versus E^/Nq under 
an AWGN environment for a smart antenna configured in 
accordance with one embodiment of the present invention. 

FIG. 9 shows simulation BER results versus Ea/N 0 under 
a Rayleigh fading environment for a smart antenna config- 
ured in accordance with one embodiment of the present 
invention. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENT 

In the following detailed description, numerous details, 
such as numbers and orientation of antennas in antenna 
arrays, are set forth in order to provide a more thorough 
understanding of the present invention. These details should 
not be understood to limit the invention. 

The present invention is believed to be particularly useful 
in CDMA (code division multiple access) systems and thus 
will be discussed in that context below. However, the 
invention is applicable to other spread spectrum systems and 
is therefore not limited to CDMA systems. The following 
description discusses a linear antenna array applied to a base 
station, however, the invention is equally applicable to the 
mobile side. The present invention is applicable to any smart 
antenna system regardless of the shape of the antenna array. 

FIG. 1 is an example receiver structure for a CDMA 
communication system with multiple user accesses. A linear 
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antenna array with M omnidirectional elements is consid- 
ered in this example. The transmitted signal from user k, 
using the baseband notation, can be written as 



s k it) = V^TottfWO exp(/(2jr/ c r + ift k )) 



(1) 5 



where P* 1 , (j)^ f c , a^t), and b^t) are the transmitted power, 
phase, carrier frequency, PN chip sequence, and information 
data bit sequence for user k, respectively. The phase ty k is 
slowly changing, compared to bit b^t), and time variable for 
ty k is dropped without loss of generality. Also, the snapshot 
rate for weight vector update and the change rate of the DO A 
angle 0 are slower than the bit and chip rate. Therefore, time 
indices for the snapshot, weight vector, and DO A will be 
also dropped. The PN chip sequence can be written as 



(2) 



20 



where a/ is the PN sequence for the k** user and P ro (t) is a 
unit rectangular pulse of chip time interval T c written as 



j 0 otherwise 



(3) 25 



The information data signal can be written as 



(4) 



x m {t) = y y aiWfri(r)or,cxi^-^^^sinffijj + 



(8) 



10 



and the complex baseband AWGN T) m (t) is given by 



15 



(9) 



Each antenna element output r m (t) after carrier down con- 
version (102) in FIG. 1 is multiplied by a complex weight 
coefficient w m , m=0,l, . . . ,M-1 (103) and then summed to 
produce the antenna array output y(t) (104), by assuming 
that an ideal low pass filter is applied to remove the double 
frequency term. Zero phase distortions at the mixers are 
assumed for simplicity. Analysis including the mixer phase 
distortions will be presented at the end of the present 
invention description. The antenna array output can be 
written as 



(10) 



30 



and 



where b * is the 1-th data bit for user k and P rt (t) is a unit 
rectangular pulse of bit time interval T b written as 



35 !Hwo, w lt ..., w M ^] T (11) 

2Hxo,x u ... f x„_J (12) 



(1 0srs7j ^ 
0 otherwise where H and T are the conjugate transpose (i.e., Hermitian) 

40 and transpose, respectively. Assume perfect PN code syn- 
chronization for simplicity. When a weight vector is 
Let k-1 denote the desired user signal index and all radio optimized, the antenna array will yield a beam pattern 
frequency mixers have the same characteristics at the towards the desired user signal and suppress interfering 
receiver. Then, the complex signal received at the m* 45 signals of different arrival angles. 

antenna element, m-0,1, . . . , M-l, can be written as The weight vector of the linear antenna array is normal- 

ized to have a unit norm, i.e., W^W^l. When the antenna 
r m (t) = V2?Tfli(/)ftj (Oaiexpf/fcr/cr + <pi - ^^anfc)) + (6) array weight vector is perfect for the desired user signal, the 



k-2 

exp^2*/ c r * <p k + 6 k - ^j^sin 9 k jj + n m (r) 



50 



weight vector can be written as 

K = 4=\l,e^ e-^-if. (13) 



55 



The antenna array response vector a is the same as the 

where m-0 is the reference antenna element (101) in FIG. 1, perfect weight vector in equation (13). Because the antenna 

d is the distance between the two adjacent antenna elements, elements are in a linear array with space d, the antenna array 

X is the wavelength of the carrier frequency, a fo 8^ and % k output y(t) would be 
are the amplitude, phase, and delay for the k fA user, 

respectively, and the bandpass noise n m (t) at the m th element 60 i ( 14) 

can be written as y(r) - XL H X.{t) = J MOM')*! +f(r) + rW 

Since coherent demodulation is assumed for the desired user 

1, Sj is set to zero. 65 where £(t) and y(t) are the interference signals due to the 

The down converted signal at the m'* antenna element is other users and the noise component, respectively, and can 

denoted by x m (t) and can be written as be written as 



06/23/2004, EAST Version: 1.4.1 



5 



US 6,483,459 Bl 



6 



M-i r k i 

m=0 I t=2 

exp^i + 5 k +■ (sin 9 X - sin O k )j j 



(15) 



(21) 



MO . 



(16) 



Then the antenna array output signal y(t) (104) can be 
written as 



10 



It is observed from equation (14) that when the weight 
vector is perfectly matched with the array response vector in 
the right hand side of equation (13), the receiver with a linear 
antenna array will have an amplitude gain of M 1/2 , i.e., the 
received signal power with M antenna elements will be M 
times stronger than that with a single antenna receiver. For 
other kinds of antenna arrays, e.g., circular or square arrays, 
etc., the same maximum antenna gain can be achieved when 
the weight vector is perfectly matched with the array 
response vector. The present invention can be applied to the 
other antenna arrays without loss of generality. The smart 
antenna output y(t) is fed to the successive post processor 
(105), e.g., demodulator, to retrieve the user data bit infor- 
mation. 

We observe that y(t) is a Gaussian random process and the 
mean of y(t) is 



m=0 

M-I K I 

( / , Indm , % 2ndm . „ X\ 
o^exd jhpt + 6 k + — — sinflj — anMJ 



(22) 



(23) 



(24) 



(25) 



where S(t), MA(t), and N(t) represent the signal, multiple 
30 access interference, and thermal noise component, respec- 
tively. The output y(t) is despread with the desired user's PN 
sequence a.,(t). The signal component of the despread output 
sampled for the bit can be written as 



35 



and the variance of y(i) is given as 



-r 

JtT b 



(26) 



/ 2ndm \ 
{n m ,M + jih»A0)txp\j — — sintf* I 



(19) 



£(«« w ,(0+/'w('))e 



4^^)]} 



40 



(27) 



t (0) 



(20) 



45 



where 



where a„ 2 -var(n m c (t))-var(n m c (t))-N 0 /2 is the variance of 
the AWGN thermal noise. 

Most of the existing smart antenna algorithms update the 50 
weight vector W=[w 0 , w,, . . . , w^J 7 " so that the weight 
vector minimizes or maximizes a cost or criteria function, 
e.g., minimizes a mean-square-error (MSE), minimizes a bit 
error rate (BER), or maximizes SINR. There are M unknown 
variables in the weight vector Instead of trying to find the 55 
M unknown variables, the present invention updates the 
DOA 0^ for user 1 and simply uses the estimate 6^ in 
equation (13) to calculate a weight vector. Only one variable 
9 j is updated instead of M variables. The present invention 
presents a method to find the optimum estimate Q 1 adap- 60 
tively to maximize correlation output with the tentative 
decision feedback bit. Two different step sizes A0 will be 
employed for the updating in the preferred embodiment. For 
the initial search, a large step size will be employed, and for 
the tracking a fine step size will be employed. Thus, the 65 
weight vector (107) using an estimated DOA 6j for user 1 
can be written as 




ai(f)flj(r)rff = n. 



When the adaptive estimation Q 1 approaches 8 lf the 
desired signal component S(lT fc ) for the X th bit can be written 
as 

The interference and thermal noise component of the 
despread output sampled for the Y h bit can be written as 
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MAtfT b ) 



= I A 



AfA(Oai(Orfr = 



— ffu<*ejcp(j(ft *-<5*)) 5^-- 

i ^pi^SrMi-"^*) 



/V(/7V) 



and 



(31) 



m-0 



where 



ffu 



=4 " 



T t )ai (r)^(r)rfr. 



(32) 



The real part of the despread output at the sampling time 
instance (l+l)T 6 is preferably used instead of the magnitude. 
Let Yl(lT fc ) denote the real part of the weighted and 
despread correlation output. Then, 



Yl{lT b yRe{S{lT b yMA{lT b yN{lT b )}. 

The MA(IT & ) and N(lT fc ) have zero means, i.e., 

E{MA(lT b )}=0 
E{N(lT b )}-0. 



(33) 

(34) 
(35) 



Let Yl denote the expectation of the real part of the 
despread output multiplied by the decision feedback bit b,\ 
Then, Yl can be written as 



Yi(B u ft^{Kl(fll)fi/} 

Y\(9 1 J l ) = E{Y/(lT b )b l l ) 

= E{Re{s{iT b )h l t + MA[lT b )h\ + N{lT b )h\}) 
= Kt{s{lT b )b\) 



1 .gi^Mi-^i) 



if 6/b/ol. Tl(B 19 6j is the taction of (sin e^sin Bj). FIGS. 
3 and 4 show the shape of Yl except the amplitude 
difference, for a beam pattern of M«6 and M=2 elements, 
respecively, and DOA«0° and DOA«60°, respectively. The 
solid line is the corresponding magnitude of E{Y1(1T 4 )6 / 1 }, 
and the dashed line is the real part of E{Yl(lT fr )b/}. The 
real part has an smaller beamwidth than that of the magni- 
tude which means that the real part has better resolution on 
the DOA. After an intitial estimation of the DOA from the 
desired user is found, the DOA tracking algorithm such as in 
FIG. 2 is activated to track the desired user. As a result, the 
beam will focus on the desired user to keep the maximum 
correlation output using the antenna output. 



If a pilot channel (106) is available in a CDMA system, 
(30) the pilot channel signal may be used as the desired signal 
may be used to track the DOA of the desired user. The pilot 
channel signal amplitude may flunctuate under fading envi- 
5 ronment. This fading environment effect will be discussed 
below. If a pilot symbol pattern is available, the pilot pattern 
b, 1 , which is known to both the transmitter and the receiver, 
may be used; or the temporary decision feedback informa- 
tion bit 6/ may be used for b, 1 . The decision feedback 
10 information bit b, 1 is the postprocessor (105) output, which 
is the binary demodulation output in the present invention. 

FIG. 2 shows an example flow chart of the proposed 
adaptive DOA tracking algorithm which is used at the smart 
antenna weight vector update (107) in FIG. 1. Here the 
15 assumption is that the angle is changing slowly, compared to 
the bit rate. This is a reasonable assumption. For example, 
suppose a mobile station is 1 km away from a base station, 
its moving speed is 100 km/hour, the data rate is 9600 bits 
per second (bps), and the DOA changes 1.5 degrees. Then 
20 (15-Jt403-3600-9600)/(180-100-103)=9047 bits will have 
been transmitted during the 1.5 degree change of DOA. For 
higher data rate, more data bits can be transmitted. Suppose 
the weight vector is updated every 20 bits. Then, 450 
updates can be made for 1.5 degree DOA change. We also 
25 assume zero angle spread for each multipath for simplicity. 
If the DOA of the desired signal is known, then the weight 
vector can be automatically obtained by using, for example, 
equation (13). However, the DOA of the desired signal is 
usually changing due to the mobile movement. Hence an 
30 adaptive DOA tracking algorithm is desirable in practice. 
Assume a power control is used so that the received power 
from any user is approximately equal to the others, which is 
a reasonable assumption in practice. The antenna array 
output will have a peak output when the beam is pointing to 
35 the desired user and the PN spread processing gain is larger 
than the number of users. 
The estimation 0! is updated as 

e^i^OyAe (4i) 

where i denotes the update index. It is not nessessary to track 
the DOA at chip or bit rate since the DOA changes slowly, 
compared to the bit rate. When B^B,, the weight vector is 
steered pointing to the desired user. 

Preferred embodiments of the initial DOA estimation 
45 process (also referred to as the DOA acquisition process) 
and the DOA tracking process will now be discussed in 
further detail. As discussed above, the initial DOA estima- 
tion process is necessary when the base station does not 
know the DOA angle for a desired user. The goal of the DOA 
50 estimation process is to estimate the DOA such that the 
estimate falls within the main lobe of the transmitted signal. 
Generally, the process comprises the following steps: form 
an initial estimate of the DOA angle, calculate a weight 
vector using this initial estimate, obtain the correlator output 
55 over some period using the weight vector, and repeat the 
previous steps until the estimated DOA is acceptable (i.e., 
within the main lobe of the transmitted signal). 

In one highly preferred embodiment, the estimated DOA 
angle is initialized to an angle corresponding to the sector 
60 boundary and then stepped (incremented or decremented, 
depending upon the initial estimate) by a coarse step A8 until 
the entire sector has been traversed. It should be noted that 
what is meant by "an angle corresponding to a sector 
boundary" includes both an angle equal to the sector bound- 
65 ary as well as an angle approximately equal to the sector 
boundary, such as an angle that is within some fraction (e.g., 
one halt) of a coarse step A0 of the sector boundary. The 
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coarse step increment A8 is preferably set to one half of the 
beamwidth of the main lobe, as shown in the fifth column of 
Table 1 below. This is done to ensure that at least one of the 
initial estimates is within the main lobe such that the 
correlator output will exceed the correlator output for other 
angles corresponding to sidelobes. The initial estimate is 
used to obtain a trial set of weight vectors. The correlator 
output is then obtained using the trial weight vectors over 
one bit (symbol) period. The estimated DO A angle is then 
updated, preferably every bit time interval, with a coarse 
step A0. This process continues until the entire sector has 
been traversed. The angle with the largest correlator output 
is then chosen as the acquired DO A and the tracking 
algorithm begins. The number of steps to search 180 degrees 
(it is not necessary to search through 180 degrees where a 
sector size is smaller, such as 120 degrees) is shown in Table 
1 for different numbers of antenna elements. At most only 
about 25 iterations (25 bit time intervals) are required to 
bring the beam pattern pointing to the desired user from any 
worst initial angle. This method is especially effective in 
systems in which power control is not used. 

Alternatively, rather than comparing correlator outputs for 
several DOA angles, comparing them to each other to 
choose the corresponding DOA angle, correlator outputs for 
estimated angles could be compared to a threshold. In such 
a scheme, the first estimated DOA angle with a correlator 
output exceeding the threshold is selected. For example, a 
first estimate DOA is made (this could be done using a 
random process, or the first estimate could correspond to a 
sector boundary as discussed above) and the corresponding 
correlator compared to the threshold. If the threshold is 
exceeded, that angle forms the initial estimate and the 
tracking process begins. Otherwise, the estimated angle is 
updated again (either randomly or by the coarse step size) 
and the process continues until the threshold is crossed. 

Table 1 shows example antenna array beam pattern 
parameters, for a 180° search range, where the number of 
antenna elements varies from 2 to 10. 

TABLE 1 

Antenna array beam pattern parameters 



Amplitude 
Number Angle Relative 

of error to the 
antenna (in maximum 
elements degrees) gain 



Number of 
steps in 
the 180° 
search range 
during the initial 
search period 



A0 (in degree) 
4 main 
lobe beamwidth x x h 
during the initial 
search period 



2 


11.0 


91.3% 


3 


60 


3 


6.0 


91.3% 


6 


31.8 


4 


4.4 


90.2% 


9 


21.6 


5 


3.3 


90.5% 


11 


16.4 


6 


2.7 


90.3% 


14 


13.2 


7 


2.3 


90% 


17 


11 


8 


1.9 


90.8% 


19 


9.6 


9 


1.7 


90.5% 


22 


8.4 


10 


1.5 


90.7% 


25 


7.4 



45 



50 



55 



When an acceptable estimated DOA is acquired, the DOA 
tracking algorithm is activated with a fine step increment AG. 
The estimated DOA during the initial acquisition period 
used as the initial value for the fine tracking period. The 
choice of fine step increment A8 depends on the number of 
antenna elements. Generally, the larger the number of 
antenna array elements is, the smaller the A0 will be 
desirable. It is a function of the number of array elements M. 
In some embodiments, the fine step increment A9 is set to 2° 
because the bit error rate is affected insignificantly by the 
beam direction error of 2° when an antenna array of M-2 



10 



elements is used. In other embodiments, the fine step incre- 
ment A8 is based on the angle error. In one highly preferred 
embodiment, the fine step increment is set to 10% of the 
angle error from Table 1 corresponding to the number of 
elements in the array. 

Reference is made to FIG. 2. The estimated DOA ange is 
initialized for two trial angles at Step 200. The following is 
a preferred updating (tracking) procedure for Q 1 at the i th 



60 



65 



Step 201: Generate two sets of weight vectors Wl and W2 
as shown in FIG. 2 with the estimated B a . Here the 
subscripts 1 and 2 represent two trial angles; early angle 
8.J-A8/2 and and late angle 6.,+A6/2, respectively. 



15 wl m = ^Le-^^i-T), »-0.1 A#-l 

w2 m = _^"1-HVt), m = 0, 1 M - 1. 



(42) 
(43) 



20 

Step 202: Measure two despread outputs Yl(lT b ) and 
Y2(lT fc ) in equation (33) at each sampling time instance 
(l+l)T b by using Wl and W2 respectively, and X in 
equation (12) and PN sequence. 
25 Step 203: Take the average of Yl(lT b )rV and Y2(lT b )b/ 1 in 
equation (36) over L bit time intervals, 



30 



1 ^ 

1 L ~ l 

Y2 = l Y 1 Y2{lT b )b l . 



(44) 



(45) 



35 Step 4: Update the estimated angle 8 j as the following 

+ if Y1<Y2 W 



di(i+0 = 



40 



8i(0, otherwise 



In step (203) the time average is used to approximate the 
ensemble average in equation (36). At step 205, the weight 
S£ is updated using the new estimated angle. 

The above-described tracking process (steps 201 to 204) 
ideally continues for as long as the incoming signal is 
received. However, under some circumstances, tracking 
may be lost, such as when the DOA may suddenly change 
due to a change in the environment, or the initial DOA 
estimate may correspond to a sidelobe rather than the main 
lobe. It is desirable to detect such events and to re-acquire 
the DOA by triggering the acquisition process. One indica- 
tion that could be used to trigger the DOA acquisition 
process is that the received correlation output after despread 
is less than a threshold, e.g., 50% of the peak gain of the 
antenna array. Another scheme is to trigger the acquisition 
process when the bit error rate or the frame error rate (or 
some other signal quality measurement) crosses a threshold. 
A third scheme is to keep the acquisition process running in 
parallel and reinitialize the tracking process if the correlator 
output for any estimated angle exceeds the correlator output 
for the current angle determined by the tracking process. 

In practical channel environments, the desired user signal 
may come from different DO As with time delays. If the 
desired signals are from the same DOA with different time 
delays, then a RAKE receiver structure can be employed to 
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/-o u i-o 



-r 



2M 



by using eqution (26). 

When the beam pattern is pointing to the DOA of the 
desired user, the desired signal output of the antenna array 
can be written as 



' m-0 



and 



7i = f{Re[5(/r fr )b;]} 



(50) 

(51) 

(52) 
(53) 



m=0 v 

. Pi sinfl 



The degradation of the gain due to the mixer phase 
distortions is about sin Q/Q. For example, when Q=60°, the 
degradation is only 10 log(0.827)— 0.825 dB. The angle 
error due to the mixer phase distortions depends on the 
number of the antenna array elements. 



12 



collect the resolvable signals. If the desired signals are from 
different DOAs, several duplicate weight vector generators 
(107 in FIG. 1) can be employed to collect the desired user 
signals received from different resovable DOAs. 
Fading 

The proposed DOA tracking algorithm discussed above is 
still valid under Rayleigh fading environment. The reason is 
that the fading amplitude of the desired user only affects the 
amplitude of the despread output at each sampling instance. 
For a given DOA, if Yl(lT b )<Y2(lT fc ), then a^YlOT^a^ 
/Y2(lT fc ), since ,>0. Even with a u changing during L bit 
time interval, the following equality still holds: 



Letip denote the equivalent DOA tracking error due to the 
mixer phase distortions, and 6/ and 6, be the true and the 
estimated DOA, respectively. Then, & 1 «0 1 +i1j. The degraded 
output Yl due to the equivalent DOA tracking error can be 
written as 



(54) 



10 



(47) 
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The equivalent DOA error due to the mixer phase distor- 
tions can be determined by the following equation: 



Hence the DOA tracking procedure remains unchanged. 
No Phase Calibrations 

In a practical system, the characteristics of the mixers are 20 
different and they introduce different phase distortions at the 
mixers in general. Let ui m be the phase distortion due to the 
m th mixer. Without loss of generality, u> OT is assumed to be 
uniformly distributed between +Q 20 ] and {\\) m } are 
independent. Including the mixer phase distortions, the 
signal component of the antenna array output after PN 
despreading can be written as 



Re 



1 _ gfeTT-bWl+VH 



- e > ? 5r w «(*!**>-«M*|. 



Afsinfl 



(55) 



25 



where op is the DOA error due to the mixer phase distortions. 
As a special case, let 9 a «0°. Then, the equivalent DOA 
tracking errors due to the mixer phase distortions are listed 
in Table 2 for M=2, 3, 5, 7 antenna elements. 

TABLE 2 



(48) 30 


Equivalent DOA tracking error due to the mixer phase distortions with 




mixer phase 


difference of ranee fi - 60° 


(49) 


M 






2 


15.85° 


35 


3 


8.6° 




5 


4.5 




7 


3.07° 



40 The tracking algorithm will work if the initial DOA guess 
is within the main lobe. The main lobe width of the beam 
pattern will affect the coarse acquisition time. The wider the 
main lobe width, the shorter the acquisition time. When the 

4 5 DOA is larger than 60° or less than -60°, the initial DOA 
search might take a wrong direction due to the poor reso- 
lution of the antenna array beam patern using the DOA. If 
the true DOA of the desired user is about 60°, a correlation 
output with a DOA guess in [-60°,-70°] might be higher 

50 than that with a guess of around 60° due to noise and MAI. 
This can be observed from FIG. 4. Fortunately a typical base 
station employs three sectors. Then, the DOAs for each 
sector can be reasonably restricted to within [-60°, 60°]. 



55 FIG. 5 shows the probabilities of correct direction detec- 
tion at the initial search versus the bit-energy-to-noise power 
spectral density ratio E^/Nq in dB for antenna arrays of 2, 4, 
and 6 elements. The N 0 represents thermal noise power 

60 spectral density in watts per Hertz. For simulation, binary 
phase shift keying (BPSK) modulation is employed. An 
antenna array of M«2 elements is used. A total of nine users 
are assumed by way of example in the system. The data rate 
is 9600 bps. The carrier frequency f c is 1.9 Giga Hertz. A 

65 Gold sequence of length 31 is used as the PN sequence. 
Under Rayleigh fading environment, the mobile speed V is 
set to 50 km/hour. 
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selecting the estimated direction of arrival angle as a 

TABLE 3 direction of arrival angle based upon the despread 

— output. 

List of the number of computations per 2. The method of claim 1, wherein the selecting step is 

snapshot for different smart antennas $ performed hy comparing the despread output to a threshold. 

Number of complex 3 ' ^ method of daim 2 > tatihst comprising the steps of: 

computations/snapshot updating the estimated direction of arrival angle if the 

~~~ estimated direction of arrival angle is not selected in the 

Adaptive algorithms Multiplication Addition selecting Step; and 

Least-Mean-Square (lms) ^ 2M + 1 2M io repeating the updating, calculating, weighting, and 

Normalized lms 3M + l 2M despreading steps if the estimated direction of arrival 

Recursive Least-Square (rls)* 2M 2 + 3M 2M 2 + 2M - l angle is not selected in the selecting step. 

Optimum adaptive convergence 4M + l 2M - l 4, The method of claim 3, wherein the estimated direction 

Mate output power criteria without 3M + 1 3M - 1 ° f 5™ 1 » by / StC u P ™°- tU . . 

Lagrange multiplier 15 5 - ^ method of claim 4, wherein the step size is 

Present invention" 2M 2M - 2 approximately one half of a beamwidth of the incoming 

d. shim and s Choi 5M + 1 3M - 1 spread spectrum signal. 

blind adaptive algorithm 5 met hod of claim 3, wherein the estimated direction 

: " ~~ z TZ " T~ of arrival angle is modified by randomly selecting a new 

*Durect RLS requires about 2M + 5M operations per snapshot. ... » • j? • , « . 

-The present invention does not update DOA every chip or bit. on diction of arrival angle within a sector corresponding to 

20 the smart antenna array. 

Table 3 compares the list of the number of computations 7. The method of claim 1, wherein the estimated direction 

per snapshot for different smart antenna algorithms where M of arrival angle is initialized to correspond to a boundary of 

is the number of antenna array elements. Table 3 shows that a sector in the initializing step. 

the present invention has a smaller number of computations 8. The method of claim 7, further comprising the steps of 

per snapshot than most of the existing smart antenna algo- ™ updating the estimated direction of arrival angle by a step 

rithms. an g le and repeating the updating, calculating, weighting, 

By updating the DOA every L bits where L-20 bits were an( J despreading, steps until the estimated direction of 

used throughout the invention, FIGS. 6 and 7 simulation test arrival angle has traversed the sector, 

results for show the transient as well as steady state of the 9 The method of claim 8, wherein the step angle is 

DOA tracking behavior under AWGN and RayLeigh fading 30 approximately one half of a beamwidth of the incoming 

environments, respectively, when M=2, Ej/N 0 -6 dB, true spread spectrum signal. 

DOA=30° and initial DOA estimation was 50°. In FIG. 7, 10- The method of claim 8, wherein the selecting step is 

mobile speed V-50 km/hour, and carrier frequency f c =1.9 performed by selecting the estimated direction of arrival 

Giga Hertz. We observe that for both AWGN and fading angle corresponding to the largest despread output from 

environments, the proposed algorithm takes only about 10 35 among all estimated direction of arrival angles in the sector, 

and 18 iterations when the mixer phase distortions are absent 11- The method of claim 1, wherein the array is linear, 

and present, respectively. 12. The method of claim 11, wherein each element of the 

FIGS. 8 and 9 show simulation BER results versus EJNq weight vector is calculated according to an equation of the 

under AWGN and Rayleigh fading environments, f° rra 

respectively, by employing the present DOA tracking 40 

algorithm, maximum output criteria without Lagrange mul- ^ _ Ce -j-^— 
tiplier (3), and Shim and Choi's blind adaptive algorithm 
(4). The simulation results indicate that the smart antennas 

configured in accordance with one embodiment of the where W^ is a weight for Mth element of the array, C is a 

present invention have the same performance, compared to 45 constant, M is an antenna array element number, X is a 

the other smart antenna algorithms. In these simulations, wavelength of the incoming spread spectrum signal, d is an 

BPSK modulation is applied. M«2 antenna elements array is antenna array element spacing, and Q 1 is the estimated 

used. A total nine users are assumed in the system. The data direction of arrival angle. 

rate is 9600 bps. The carrier frequency f c is 1.9 Giga Hertz. 13. The method of claim 11, wherein the weight vector is 

Walsh sequences of length 32 are used as the PN sequences. 50 calculated according to an equation of the form 
Under Rayleigh fading environment, the mobile speed V is 
set to 50 km/hour. 
What is claimed is: 

1. A method for acquiring a direction of arrival angle for 

an incoming spread spectrum signal comprising the steps of: 55 

initializing an estimated direction of arrival angle of an ^ herein ^ * th f wei 8 ht vector > M * 311 antenna arra y 

incoming spread spectrum signal; element nllmber > K 15 a wavelength of the incoming spread 

.... ... „ spectrum signal, d is an antenna array element spacing, 6, is 

calculating a weight vector for a smart antenna array ^ of arriya , J k md ^ ^ ^ 

based on the estimated direction of arrival angle, the 6Q 

weight vector being matched to an array response u The memod of daim x wherein aR im in tion 

vector for the antenna array; of ^ dcspfead Qutput fa ^ u$cd {n ^ gtep 

weighting signals received by the smart antenna array 15. The method of claim 1, wherein both a real and an 

with the weight vector to form a weighted output imaginary portion of the despread output is used in the 

without performing phase calibration; $5 selecting step, 

despreading the weighted output to form a despread 16. The method of claim 1, wherein the incoming spread 

output; and spectrum signal is a code division multiple access signal. 
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17. The method of claim 1, wherein the incoming spread 28. The method of claim 27, wherein each element of the 
spectrum signal is a pilot channel signal. first weight vector is calculated according to an equation of 

18. The method of claim 17, wherein the despread output the form 
is multiplied by a pilot symbol pattern. 

19. The method of claim 1, wherein the incoming spread 5 p***^ 
spectrum signal is multiplied by a tentative decision feed- ** B e 

back bit. 

20. The method of claim 1, wherein the signals weighted where is a weight for Mth element of the array, C is a 
in the weighting step are collected over a bit period. constant, M is an antenna array element number, X is a 

21. A method for tracking a direction of arrival angle 10 wavelength of the incoming spread spectrum signal, d is an 
comprising the steps of: antenna array element spacing, and Q x is the first trial angle. 

initializing an estimated direction of arrival angle of an 29. The method of claim 27, wherein the first weight 

incoming spread spectrum signal; vector is calculated according to an equation of the form 

forming a first trial angle, the first trial angle being spaced 

apart from the estimated direction of arrival angle by a 15 _ l r e -j^iurf, g j* < y- i >« B g i r T 

first angular difference in a first direction; ~~ ~ VaT ' ' e ' " * C -I 

calculating a first weight vector for a smart antenna array 

based on the first trial angle, the first weight vector wherein ^ is ^ first wd ht vcc M is an antenna array 

being matched to an array response vector for the 2Q e i ement numberA is a wavelength of the incoming spread 

*' spectrum signal, d is an antenna array element spacing, 8^ is 

weighting signals received by the smart antenna array the first trial angle, and T is the transpose. 

with the first weight vector to form a first weighted 30. The method of claim 21, wherein an imaginary portion 

output without performing phase calibration; 0 f the first despread output is not used in the selecting step. 

despreading the first weighted output to form a first 25 31. The method of claim 21, wherein both a real and an 

despread output; and imaginary portion of the first despread output is used in the 

selecting the first trial angle as a new direction of arrival selecting step. 

angle based upon the first despread output. 32 - The method of claim 21, wherein the incoming spread 

22. The method of claim 21, further comprising the steps spectrum signal is a code division multiple access signal, 
of: 30 33. The method of claim 21, wherein the incoming spread 

forming a second trial angle, the second trial angle being * S \^°} C *l nnt } sig ? aL 
spaced apart from the estimated direction of arrival . 34. The method of claim 33 wherein the despread output 
angle by a second angular difference in a second ^ mulched by a pilot symbol pattern 
direction, the second direction being different from the 35 * ™ c mcth ? d of daim 21, wherein the incoming spread 
first direction* 35 s P ectrum sign**! I s multiplied by a tentative decision feed- 
calculating a second weight vector for a smart antenna M tu ^ f ^ ia u - i i ■ 

, & , • 1 i.i. ■ t_ 4 36. The method of claim 24, wherem a loss of tracking is 

array based on the second trial angle, the second weight { a b ^ for 

vector being matched to an array response vector for spread ^ tQ & lhreshold 

the antenna array, 4Q 3? ^ method of claim ^ wherein a loss of tracking is 

weighting signals received by the smart antenna array detected by comparing a frame error rate for the incoming 

with the second weight vector to form a second spread spe ctrum signal to a threshold, 

weighted output; and 38 ^ method of claim 2 ^ wherein a loss of tracking is 

despreading the second weighted output to form a second detected by comparing a signal strength to a threshold, 

despread output; 45 39. A method for acquiring a direction of arrival angle for 

wherein the selecting step is performed by selecting the an incoming spread spectrum signal comprising the steps of: 
first trial angle as a new estimated direction of arrival initializing an estimated direction of arrival angle of an 
angle if the first despread output is larger than the incoming code division multiple access spread spec- 
second despread output, selecting the second trial angle t rum signal to correspond to a boundary of a sector; 
as a new estimated direction of arrival angle if the so calculati a wei ht vector for a smart 
second despread output is larger than the first despread based Q * ^ of ^ x } ^ 
outpu , and otherwise leaving the estimated direction of we{ fat vedor be{ matched {Q aQ ^ 
arnval angle unchanged vector for the antenna array; 

23. The method of claim 21, further co mpnsing the steps ., . , . , 

of repeating the forming, calculating, weighting, despread- 55 wei S htin g "go* received by the smart antenna array 

ing and selecting steps until a loss of tracking is detected. over u a s y mbo1 P™*? ™& ih f ™ ctor to form a 

24. The method of claim 22, further comprising the steps weighted output without performing phase calibration; 
of repeating the forming, calculating, weighting, despread- decorrelating the weighted output to form a decorrelated 
ing and selecting steps until a loss of tracking is detected. output; 

25. The method of claim 21, wherein the selecting step is 60 multiplying the decorrelated output by a symbol bit corre- 
performed by comparing the first despread output to a sponding to a symbol represented by the incoming code 
threshold. division multiple access spread spectrum signal during the 

26. The method of claim 21, wherein the first angular symbol period to form a despread output; 

difference is approximately equal to ten percent of an angle changing the estimated direction of arrival angle by a step 

error of the smart antenna array. 65 size approximately equal to one half of a beamwidth of 

27. The method of claim 21, wherein the smart antenna the incoming code division multiple access spread 
array is linear. spectrum signal and repeating the calculating, 
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weighting, decorrelating, and multiplying steps for 
each estimated direction of arrival angle until the sector 
has been traversed; 
selecting the estimated direction of arrival angle corre- 
sponding to the largest despread output as the direction 
of arrival angle. 
40. The method of claim 39, wherein the array is linear 
and the weight vector is calculated according to an equation 
of the form 



[l.e-'T-^l, 



15 



30 



35 



wherein W is the weight vector, M is an antenna array 
element number, X is a wavelength of the incoming spread 
spectrum signal, d is an antenna array element spacing, B a is 
the estimated direction of arrival angle, and T is the trans- 
pose. 

41. The method of claim 40, wherein only a real portion 
of the decorrelated output is multiplied by the symbol bit to 20 
form the despread output. 

42. The method of claim 41, wherein the incoming spread 
spectrum signal is a pilot channel signal and the symbol bit 
corresponds to a pilot symbol pattern. 

43. The method of claim 41, wherein the symbol bit 25 
corresponds to a tentative decision feedback bit. 

44. A method for tracking a direction of arrival angle 
comprising the steps of: 

initializing an estimated direction of arrival angle of an 
incoming spread spectrum signal to a value within a 
main lobe of the incoming spread spectrum signal; 

forming a first trial angle, the first trial angle being spaced 
apart from the estimated direction of arrival angle by a 
step angle in a first direction; 

forming a second trial angle, the second trial angle being 
spaced apart from the estimated direction of arrival 
angle by the step size in a second direction; 

calculating first and second weight vectors for a smart 
antenna array based on the first and second trial angles, 
respectively, the weight vectors being matched to an 
array response vector for the antenna array; 

weighting signals received by the smart antenna array 
during a symbol period with the first and second weight 
vectors to form a first and second weighted outputs 45 
without performing phase calibration; 

decorrelating the first and second weighted outputs to 
form first and second decorrelated outputs; 

multiplying a real portion of the first and second decor- 
related outputs by a symbol bit corresponding to a 50 
symbol represented by the incoming code division 
multiple access spread spectrum signal during the sym- 
bol period to form first and second despread outputs; 

selecting the first trial angle as a new estimated direction 
of arrival angle if the first despread output is larger than 55 
the second despread output, selecting the second trial 
angle as a new estimated direction of arrival angle if the 
second despread output is larger than the first despread 
output, and otherwise leaving the estimated direction of 
arrival angle unchanged; and 



repeating the calculating, weighting, decorrelating, 
multiplying, and selected steps until a loss of tracking 
is detected. 

45. A method for tracking a direction of arrival angle for 
a signal comprising the steps of: 

estimating an estimated direction of arrival angle; 
calculating a weight vector for an antenna array compris- 
ing a plurality of antenna elements by matching the 
1 0 weight vector to an antenna array response vector based 
on the estimated direction of arrival angle; 
obtaining a signal weighted by the weight vector from the 

antenna array without performing phase calibration; 
updating the estimated direction of arrival angle based on 
a parameter of the signal. 

46. An apparatus for acquiring a direction of arrival angle 
for an incoming spread spectrum signal, the apparatus 
comprising: 

a processor connectable to a smart antenna array, the 
processor being configured to perform the steps of 
initializing an estimated direction of arrival angle of an 

incoming spread spectrum signal; 
calculating a weight vector for the smart antenna array 
based on the estimated direction of arrival angle, the 
weight vector being matched to an array response 
vector for the antenna array; weighting signals 
received by the smart antenna array with the weight 
vector to form a weighted output without performing 
phase calibration; 
despreading the weighted output to form a despread 
output; and 

selecting the estimated direction of arrival angle as a 
direction of arrival angle based upon the despread 
output. 

47. The apparatus of claim 46, further comprising a smart 
antenna array connected to the processor. 

48. An apparatus for tracking a direction of arrival angle, 
4Q the apparatus comprising: 

a processor connectable to a smart antenna array, the 
processor being configured to perform the steps of 
initializing an estimated direction of arrival angle of an 

incoming spread spectrum signal; 
forming a first trial angle, the first trial angle being 
spaced apart from the estimated direction of arrival 
angle by a first angular difference in a first direction; 
calculating a first weight vector for a smart antenna 
array based on the first trial angle, the first weight 
vector being matched to an array response vector for 
the antenna array; 
weighting signals received by the smart antenna array 
with the first weight vector to form a first weighted 
output without performing phase calibration; 
despreading the first weighted output to form a first 
despread output; and selecting the first trial angle as 
a new direction of arrival angle based upon the first 
despread output. 
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